Singular behavior of an entangled state for a one-dimensional quantum spin system 
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We studied the entangled state for a one-dimensional S = 1/2 antiferromagnetic quantum spin 
chain in a transverse field. We calculate the ground state using the density matrix renormalization 
group and discuss how the entangled state changes around a quantum phase transition (QPT) 
point. By analyzing concurrence C(p) for two-qubit density matrix p after the Lewenstein-Sanpera 
decomposition, p = Ap s + (1 — A)p e , where p s is a separable density matrix and p e is a pure 
entangled state obtained by a linear combination of Bell states, we find singular behaviors both 
in C(p e ) and 1 — A at the QPT point.® C(p e ) includes the effects of quantum fluctuations, which 
manifest the competition between the antiferromagnetic spin fluctuation and the effect of transverse 
field in the transverse Ising model. The quantum fluctuation shows the singular maximum at the 
QPT point as expected from the general picture of phase transition. In contrast, 1 — A reveals the 
singular minimum at QPT point. 



Recently, the entanglement properties of quantum spin 
systems have attracted much interest P,0,0,SIEIEHI3 
. One of the most interesting topics is the behavior of en- 
tanglement at various types of quantum phase transition 
(QPT) that are induced by the change of an external 
field or coupling constant at absolute zero temperature. 
This is because we can expect that the quantum fluctu- 
ation around QPT points triggers anomalous behavior. 
To study the entanglement embedded in the quantum 
fluctuation, it is useful to analyze the concurrence || for 
two-neighboring spins of the ground state. In practice, 
Osterloh et al. have shown that the derivative of the 
concurrence with respect to a transverse magnetic field 
diverges at the transition point [1| . However the concur- 
rence itself does not show a maximum. At a glance, this 
seems to go against the general picture of phase transi- 
tion because the fluctuation becomes the maximum at 
the critical point in general. The reason the maximum 
point shifts from the critical point is not known precisely. 

To give an answer to the above problem, in this 
per, we apply the Lewenstein-Sanpera decomposition 
to analyze the pairwise-entanglement of one-dimensional 
(ID) quantum spin systems. Lewenstcin and Sanpera 
have shown that any bipartite density matrix p can be 
represented optimally as a sum of a separable state p s 
and an entangled state p e . For two-qubit systems, they 
have also shown that the decomposition reduces to a mix- 
ture of a mixed separable state and an entangled pure 
state (p e = |* e )(* e |) as p = Kp s + (1 - A)p e , thus all 
non-separability content of the state is concentrated in 
the pure entangled state. This leads to an unambiguous 
measure of entanglement for any two qubit state as en- 
tanglement of the pure state |\E , e )(4' e | multiplied by the 
weight (1 — A) of pure part in the decomposition. 

We have employed the density matrix renormaliza- 
tion group (DMRG) method 11] to calculate the ground 
state of the Hamiltonian (|TJ of ID system. In particu- 
lar, we make use of the infinite DMRG method in or- 
der to compute bipartite density matrixes in the bulk 
limit. The infinite DMRG enables us to treat ID sys- 
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FIG. 1: The uniform and staggered magnetization: (a) Ising 
model (Jy = 0, J± = 1) and (b) XY model (Jy = 1, Jj_ = 0). 



terns with high accuracy even around QPT point |13| . 
But, for the Ising case (Jii = 0), since the exact ground 
state was obtained by previous works [12j , we have taken 
advantage of the exact solution to check our numerical 
results. Our numerical results clarify that (i) the concur- 
rence C(p) for the pure entangled state |^ r e )(^ r e | shows a 
singular maximum at the QPT point and (ii) the weight 
(1 — A) of the pure part shows a singular minimum at 
the point. This is the reason that the maximum point of 
C{p) = (1 - A)C(p e ) shifts from the QPT point. 

The Hamiltonian for the ID S = 1/2 anisotropic 
Heisenberg model in a transverse field is given by 

h = y][j\\ (Sfsf +1 + sfsf +1 ) + J_i_sfsf +1 ] 

i 

-hJ2 S ? (!) 

i 

where Sf is the 5=1/2 spin operator at site z, and 



2 



h is the strength of the transverse magnetic field. The 
exchange coupling constant in the easy-plane (xy-pl&ne) 
is denoted by Jy (> 0) and that of the easy-axis (2-axis) 
isJj_(>0). 

The properties of Hamiltonian have been studied 
for several years. When h = 0, the ground state of this 
system is described as follows. For Jj| < J±, the systems 
belong to the Ising universality class, then the property of 
the systems are characterized with the spontaneous stag- 
gered magnetization | (S z ) | that has a finite value even at 
h = 0. 

On the other hand, for Ju > J±, the systems belong 
to the Tomonaga-Luttinger liquid phase, therefore there 
is no spontaneous staggered magnetization. The prop- 
erty of this phase is that the correlation length £ of spin 
fluctuation diverges in the whole region. Namely, in ID 
5 = 1/2 spin systems, the critical behavior appears in 
the Tomonaga-Luttinger liquid phase and in the phase 
transition points as well. 

Then, it should be mentioned regarding, magnetic 
properties, there is a special value of the transverse field 
called a classical point h c i. This point is exactly known 
as h c i = ^/2 J|| ( J|| + JjJ, and the ground state can be 
expressed by a simple product state 0, ^| . 

In Fig. ^ we show the uniform and staggered mag- 
netization for the Ising model and XY model. The 
classical point corresponds to h = and h = ^/2, re- 
spectively. Here, two magnetizations have the relation 
\J{S X ) 2 + (S z ) 2 = 1/2 at these points. 

To analyze the pairwise-entanglement of two- 
neighboring spins, we prepare the two-qubit density 
matrix from the wave function of the ground state ob- 
tained by DMRG. In S = 1/2 systems, Wootters has 
shown that the pairwise-entanglement (entanglement of 
formation) of a mixed state can be defined as 



C(p) — max[0, Ai — A2 — A3 — A4 



(2) 



where the X a are the square roots of the eigenvalues of 
the matrix R = p{cr y ®a y )p* [a y ®a y ) in decreasing order. 
The concurrence varies from C — for a separable state 
to C — 1 for a maximally entangled state. 

Moreover, according to Lewenstein and Sanpera [Icj . 
any two-qubit inseparable density matrix p can be de- 
composed as 



Ap s + (1 - A)p e 



(3) 



where p s is a separable density matrix and p e = \^ e ) (^e\ 
is a pure entangled state obtained by a linear combination 
of Bell states, (|V±) = ^| Tl) ± | IT), \<f>±) = 4s I TT 

>±IU» 



I * e )=a|V- > + *#+) + <#_)+<#+) 



(4) 



where |a| 2 + 



1. 



The Lewenstein-Sanpera (L-S) decomposition of a 
given density matrix p is not unique and, in general, there 
is a continuum set of decompositions to choose from. The 
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FIG. 2: Numerical results of L-S decomposition for the trans- 
verse Ising model ( Jy = 0, J± — 1) as a function of the trans- 
verse magnetic field h. h = 0.5(= h p ) denotes the QPT point 
from the AFM phase to the PM phase. 



optimal decomposition is, however, unique for which A 
is maximal. To determine the maximum A, we perform 
calculations as follows. First, we make 



{p-(l-%}/A. 



(5) 



with variables (A, a, b, c, d). Then we make use of the 
positive partial transpose (PPT) 0, 0] to determine 
whether p s is separable (p s = '^2 i PiPA ® Pb) or not. We 
find the optimum A, which is the maximum value under 
the condition that p s is separable. 

If A is optimum, the concurrence C{p) for p is related 
to 1 — A and p e of the L-S as 



C(p) = (l-A)C(p e ) 



(6) 



where C(p e ) = \a 2 - b 2 - c 2 + d 2 \. 

In Fig. EJa) and Fig. we show the values of C(p), 

C(p e ) and 1 — A for the Ising model and XY model, re- 
spectively. For the Ising case, h p = 0.5 corresponds to 
the phase transition point. In this result, we find C(p e ) 
and 1 — A show singular behaviors at h p . Also, for the 
XY model, C(p e ) and 1 — A has the singular maximum 
and minimum at h p , respectively. As a consequence of 
these contrasting behaviors, the maximum point of con- 
currence C{p) shifts from the QPT point. 

We consider the detail of the C(p e ), which is the con- 
currence of a pure entangled state. In Fig. E^b) and 
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FIG. 3: Numerical results of L-S decomposition for the XY 
model (J|| = 1, J x = 0). h ~ 1.458(= h p ) denote the QPT 
point from the AFM phase to the PM phase and h — \/2(= 
ftci) is the classical point, which can be expressed by a simple 
product state. Note that the Tomonaga-Luttinger liquid state 
appears at h = 0. 



Fig. 0{b), w e show the ratio of Bell states. The \a\ 2 indi- 
cates the probability of the \ip-} state. This probability 
vanishes entirely at h p in both cases. Therefore, we can 
consider this quantity to be the parameter characterizing 
the QPT from the antiferromagnetic (AFM) phase to the 
paramagnetic (PM) phase. Also this is the order parame- 
ter like the staggered magnetization. On the other hand, 
|c| 2 means the effect of ferromagnetic fluctuation, and 
|c| 2 = for these cases. But, if we consider the pairwise 
entanglement of next-nearest neighbor spins, |c| 2 has fi- 
nite values and behaves itself like the order parameter, 
substituting for \a\ 2 . 

In Fig. C2b), the point of h = is the Tomonaga- 
Luttinger liquid state with zero staggered magnetization. 
In the low-field region, the system is under the influence 
of the AFM quantum fluctuation due to the XY-term 
(J\\ £[Sf Sf+i + SfSf +1 }), while, in the high-field region, 
the effect of the transverse-field h becomes dominant. In 
the latter case, the effect of transverse-field is reflected in 
\b\ 2 and |g?| 2 . Here it should be noted that the total-S 12 



is not a conservative quantity for the anisotropic models 
in the transverse field. Therefore, we can regard \b\ 2 and 
\d\ 2 to reflect the spin-number fluctuation. In particular, 
\b\ 2 increases in the middle-field region. This means that 
the enhancement of C(p e ) at h p is attributed to the spin- 
number fluctuation. 

Next we focus on the classical point /i c ;. However, we 
can not determine the set of (A, a, 6, c, d) uniquely at h c i, 
because a given density matrix p is already separable. 
Therefore, we discuss the classical point from the limit of 
both sides of h c \. As |c| 2 = and |6| 2 > |d| 2 in all regions, 
C{ Pe ) is represented by C(p e ) = \a\ 2 - (\b\ 2 - \d\ 2 ). The 
first term and second term imply the effect of the AFM 
fluctuation and the spin-number fluctuation, respectively. 
This equation means that the two quantum fluctuations 
conflict with each other. Actually, in Fig. [3J C(p e ) goes 
to zero at h d = y/2, and then \a\ 2 -► 1/4 , \b\ 2 -►1/2 
and \d\ 2 — ► 1/4, respectively, while, for the Ising model, 
the classical point is h = 0. Then, the entangled state is 
approached as 



(7) 



Here the ± indicates that the ground state approaches 
11) or | due to the spontaneous symmetric break- 
ing. 

We discuss the behavior of 1 — A at around h p . As 
shown in Fig. Ufa) and Fig. E[a), 1 — A exhibits the sin- 
gular minimum at h p . Moreover, as shown in Fig. ffi a). 
1 — A shows another minimum at the point of ft. = 0, 
where the system behaves as Tomonaga-Luttinger liq- 
uid. Note that these two points correspond to the crit- 
ical point of the system, where the correlation length £ 
of spin fluctuation is infinite. These facts suggest that 
the divergence in the number N of the correlated spins 
brings the minimum in 1 — A observed at the critical 
points. In practice, an amount of the pairwise entangle- 
ment decreases necessarily when the number of entangled 
qubits increases. 

In conclusion, we have studied the pairwise- 
entanglement in ID S = 1/2 AFM quantum spin system 
in a transverse field. We have found that not only the 
derivative of the concurrence, which have shown by Os- 
tcrloh, et. al. [l|, but also the decomposed components 
indicate the singular behaviors at QPT point. We also 
expect these components show the scaling behaviors in 
the vicinity of the phase transition point. In particular, 
the probability of \ip-) is suitable for describing QPT 
from the AFM to the PM phase, because it vanishes en- 
tirely at QPT point. 
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